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We present a variational approach to evaluate relativistic corrections of order a^ to the Bethe 

logarithm for the ground electronic state of the Coulomb two center problem. That allows to 

estimate the radiative contribution at ma^ order in molecular-like three-body systems such as 

C^ hydrogen molecular ions H^ and HD^, or antiprotonic helium atoms. While we get 10 significant 

^~~^ digits for the nonrelativistic Bethe logarithm, calculation of the relativistic corrections is much 

^—^ more involved especially for small values of bond length R. We were able to achieve a level of 3-4 

^^ significant digits starting from R — 0.2 bohr, that will allow to reach lO"'^" relative uncertainty on 

5-H transition frequencies. 
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Introduction 



Considerable effort is currently devoted to high-precision laser spectroscopy of three-body molecular (or molecule- 
like) systems such as HD+ 1, 2,, Hj ;3J and antiprotonic helium 4J. These experiments aim at improving the present 
O accuracy of the electron-to-proton and -antiproton mass ratios, for which the uncertainty of spectroscopic data, as well 

I as of theoretical calculations of transition frequencies, should reach a level of about 0.1 ppb. Systematic evaluation of 

H leading QED corrections up to the ma^ order has improved the theoretical precision in hydrogen molecular ions [5] and 

Q antiprotonic helium [^ to 0.3-0.4 ppb and 1 ppb, respectively. The main source of theoretical uncertainty is the ma^ 

order one- loop self-energy correction [Tj, which so far has been evaluated only in hydrogenlike systems. Considering 
the aimed for accuracy, it is enough to calculate the relativistic Bethe logarithm with 3-4 significant digits and thus it 
fj may be obtained in the framework of the adiabatic approximation, i.e. for an electron in the field of two fixed nuclei. 
' 'r? ^'^^ hydrogen-like ions, the one-loop self-energy contribution to the binding energy of an electron is traditionally 



^^ expressed as follows [H] 

Dh A£;i_ioop = H^i^ I U4i(n) In [{Za)-^] + A,o{n)] + (Za)A5o(n) 



(1) 



^ +{Zay Ae2in) ln'[{Za)-'] + Aei{n) In [(Za)-'] + ^goW + • • ■ ■ 

> L J J 

m It is known that among ma^ order terms, the Bethe logarithm (which appears in the low-energy part of the nonlog- 
'^ arithmic contribution A40) is the most difficult quantity for numerical evaluation. Similarly, at the ma^ order the 

^? low-energy part of Aqq contains the relativistic Bethe logarithm [5HT^ which gives rise to even more severe difficulties. 
. In the present work, we describe a numerical method which allows to obtain these quantities with very good accuracy 
^r for a two-center system. 

The paper is organized as follows. In Sec. |l) we briefly outline the origin of Bethe logarithm contributions in a 
nonrelativistic quantum electrodynamics (NRQED) approach [T31[T1] and give their precise definition. In Sec.lllj the 



. . asymptotic behavior of the integrands is derived. In Sec. Ill the numerical method is described in detail, and finally, 

^ the nonrelativistic and relativistic Bethe logarithms are calculated for the hydrogen atom, hydrogen molecular ions 
and antiprotonic helium. 



I. BETHE LOGARITHM: DEFINITIONS 

In this Section natural relativistic units {h = c = m =1) are used, while starting from Sec. Ill] we switch over to the 
atomic units, which are more suitable for numerical calculations. 
As a starting point, we take the nonrelativistic Hamiltonian 

H=f+V, F = -^-^ (2) 

2m ri r2 

where ri and r2 are the distances from the electron to nuclei 1 and 2, respectively. The case of Zi = Z2 = 1 corresponds 
to the hydrogen molecular ions and Zi=2, ^2 = — 1 to the antiprotonic helium atom. 



A. The NRQED one-loop self-energy at the ma^ order (low photon energy) 

The leading order NRQED interaction with the magnetic field is determined by 

m Zm 

The first term in this expression is the "dipole" interaction, while the second one is called Fermi's interaction. 



FIG. 1: The NRQED diagram for the leading order one-loop self-energy contribution 

R may be shown [15] that the NRQED diagram in Fig. fl] with the Fermi- type interactions on one or both sides of 
the transverse photon line give vanishing contributions. Thus the low energy part contribution, which stems from the 
NRQED diagram in Fig. [T] with two dipole vertices may be expressed: 

^'- = I^ L. ? {'" - ^) (=■""' ( J^^^ki) <"■-'") - '- <*!*"> <^' 

where Sm is a " mass renormalization" term. Here and throughout this paper it is assumed that in (...) on the left 
and right-hand sides of brackets stands ipo, a stationary Schrodinger eigenstate of the Hamiltonian operator of Eq. (l2|, 
if not otherwise stated. 

To get the leading ma^ order contribution, one neglects retardation, replacing the exponential factors in (l3|) by 
unity, which leads to the nonrelativistic dipole approximation 

with 

Pnd{k) = (p{Eo-H-k)-'py (5) 

The integral in Q contains a linearly divergent term which corresponds to the electron's mass renormalization, as 
was shown by Bethe in 1947 [T^. R also contains a logarithmic term, where the dependance on the cutoff parameter 
A is canceled by the logarithmic contribution from the high-energy part P71. After these two terms are dropped out, 
the remaining nonlogarithmic contribution at order ma^ may be written [Eh = ma^ is the Hartree energy) 



Jo \ \Eo-H-k k)^l Je^ k yEo-H-k / 

and determines the numerator of the Bethc logarithm, while the denominator is expressed by 

Vin;R) = l{AV) (7) 

The Bethe logarithm itself is defined as the ratio 

Pnr{n;R)^^. (8) 

Here n denotes a set of state quantum numbers. 

B. One- loop self-energy contributions at the ma^ order 

Here, for convenience of reading, we keep the notation of \TT wherever possible. There are three types of relativistic 
corrections to the leading order expression Q, which give a contribution at order ma^ : 




FIG. 2: NRQED diagrams for relativistic corrections to the electron line, which contribute to self-energy at order ma^ . 
1. Relativistic corrections due to the Breit-Pauli interaction (Fig. [2]) 

where 



fp (Eo-H- k)-' (Hb - (Hb)) (Eo-H- k)-^ p) 



Here Q is a projection operator: Q = /— |?/'o)(^o|- Eq. (10 1 represents the third order term in the Rayleigh-Schrodinger 



perturbation theory. The relativistic Breit-Pauli Hamiltonian for the two-center problem is expressed 

Hb = -^ + ^ [4^^i5(ri) + 47TZ2S{r^)] , 

the spin interaction is neglected. 

For reasons which will be discussed later, it is convenient to split Pre into two parts: 

P^l^\k) = (^p{Eo-H- k)-' (Hb - {Hb)) {Eo - H - ky' p) (11a) 

P^l'\k) = 2(^HBQ{Eo-H)-^Qp{Eo-H-k)-'p) . (lib) 

2. Modification of vertex interactions in the self-energy diagram 





FIG. 3: The NRQED diagrams for the self-energy with modified vertex interactions at order ma' 
The next order NRQED interactions {^a^Hj ) with a magnetic field are determined by [CTfTH] 
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thus modifying vertex functions as shown on the diagrams in Fig. ([3]); one gets 

where S,P = —p^p^ — ^a^^VW, and 

P^? (k) = ( (-p'f - la^'^'v] {Eo-H- k)-' pA . (13) 



3. It remains to consider the effect of retardation (see Eq. ([3|)). We obtain the nonrelativistic quadrupole contri- 
bution, which results from the Taylor series expansion of e**^'""^ = l-l-i(kr) — (kr)^/2 + . . . 

2a^ r^ 
P^L2^^—^j kdkP,,q{k) (14) 






-(fin-rf{Eo-H-k)-'p' 



(15) 



where k — fcn. 

Similarly to the nonrelativistic Bethe logarithm considered above, the relativistic Bethe logarithm corresponds to 
the finite part of the integrals ^, ( 12 ) and ( 14), i.e. divergent terms in A must be subtracted jljj. It is thus essential 
to study the asymptotic behavior of the integrands in the fc — >■ cx) limit. 



II. ASYMPTOTIC BEHAVIOUR OF THE INTEGRANDS AT fc ^ cx) 

Our approach to obtain asymptotic expansions stems from ideas first formulated by C. Schwartz in [19J. The first 
step is to note that in expression Q the integrand's form is that of a second-order perturbation correction. It may 
be calculated via the first-order perturbation wave function -01, which can be obtained by solving the differential 
equation 



and then one calculates the integrand by evaluating 

Pnd{k) = (0o|V|Vl) . 



(16) 



(17) 



A. Nonrelativistic Bethe logarithm 



The first order nonrelativistic perturbation wave function for A; — > oo to a good extent may be approximated (see 
[ISIEO]) by 



V'i(r) 



Ziri , Z2r2 



ri 



r2 



Mr) - ^ (^ [1-e-^'-^ (l+/xri)] + ^ [l-g-^''^ (l+^^r^)] \Mr), (18) 



here fi = v 2fc. Th is fu nction has a proper behaviour at r ^ 0. 

By substituting (18) into expression j: (V^) — ]: {ipo\[H,'V]\ipi) (see for details p ^ l20 j ) one gets for the nonrelativistic 
dipole term 



1 



1 



P„.(fc) = r(V^) + ;^(AF) 



2fc2 



-^{[zy2fc-Z3lnfc] 47r(<5(ri)) 



ZiV2k-Zl\nk 



in 



Wr2))} 



(19) 



in which all terms (including the last one in Infc/fc^) are correct. For higher-order terms in 1/fc we use the same 
expansion as for the hydrogen atom 



yQji^ + Q^fr>k + Q[ 



3m 



hm-{-3 



(20) 



which is in the latter case known analytically [H] (See also [5]). 



B. Relativistic Bethe logarithm 



Substituting again the wave function tpi from (18) into the matrix elements which appear in the integrands Pi{k) 
in Eqs. (15), (11a), (lib) and (13l, one gets 



Pn,{k)^-l{^')-l 



(V^) 



+ 27r{Z,{6iv,))+Z2{Sir2))) 



2ZHV2k -Z.lnk)^^ ^^^^^^^ ^ 2ZUV2k Z^lnk) ^^ ^^^^^^^ 



F„. 



A R 

-^^nq -'-'nq 



fc2 
fc3/2 ^ /,2 



L-nq In k J-^nq 



(21a) 



P^c^'ik) . -^4. (^(r,)) - ^ 4. (^(r,)) 



4fc2 



4fc2 



R(la) 

fc3/2 



z? 



(la) 



fc2 



(21b) 



^^"^(fc) = t: {(HB-iHs)) (Eo-Hr'V') 



Zf{2V2k + Zilnk) 



A 



(lb) 



r(16) 
Ore 

J3J2 



4/c2 

fc2 



D 



47r(<5(ri)) 

(16) 



Z^{2V2k + Z2\nk) 
4fc2 



ATT{d{r2)) 



(21c) 



pyjik) 



(p^) Zf{-VSk + Zilnk) 



An ((5(ri)) 



A 



(2) 



B. 



(2) 



O/'c in "J j-^rc 



1,3/2 



B 



B 



Zi{-VSk + Z^lnk) 



47r(<5(r2)) 



(21d) 



For higher order terms, in the case of Pnq, Pre"' and Pre the form of the asymptotic expansion is found to be 



similar to the nonrelativistic Bethe logarithm (Eq. (20)), for example: 

j(la) T^ila) 



P^l"Hk) 



Br 



D. 



1,3/2 



fc2 



771 — 1 



la) 



^m+2 



(22) 



5(2) 



with equivalent expressions for Pnq{k) and Pre (k). The Pre term has an essentially different asymptotic behavior: 



(ifc) 



pZ'Kk) 



Alb) 

-TLrC 



B 



(16) 



Ot^c -^n '^ -L-^TC 






/c3/2 



fc2 



/c2 



fc 



^EE 



7n— 1 n— 



J^m/2 



(23) 



This is one of the reasons why the Pre term has been separated in two contributions. 

In actual calculations coefficients of the asymptotic expansion A, B, C, and F are calculated from expectation 
values of the operators appearing in the Eqs. (21i)-(21i), while the unknown coefficients D, Q, and S are obtained 
by fitting of the numerically evaluated integrand using Eqs. (22|-(23). 



C. Final expression of the relativistic Bethe Logarithm 



In view of the asymptotic expansion obtained in the previous paragraph, the relativistic Bethe logarithm, which is 
given by the finite part of the integrals (l9])-(14), can be written as follows: 



£ = p^^^ + /?f ) + /32 + /33 , 



(24a) 



o(a) _ 2 
^1 "3 



k dk 



py^^Hk) 



B 



(la) 



k dk 



fc3/2 



Urc 
/c3/2 






(24b) 



/3}' 



(&) 



3 7o 



kdk 



Hl'Hk) 



A 



(lb) 



B 



(lb)' 



k dk 



k fc3/2 

P(f)(fc)- 



J-i-rc 



k 



/fc3/2 fc2 p~ 



(24c) 



/32 



3 7o 



kdk 



A R 



fc dk 



k fc3/2 

4 R 

p r/c^ _ p - ^^ - =^ 



^nq in "- ^ng 



(24d) 



/33 = 



kdk 



p^m- 



4(2) 0(2) 

~k p72 



fc dk 



PlVik) 



A. 



(2) 



W^ fc2 ^ 



(24e) 



Here the terms which are subtracted from the first line of each equation appear in the expansions due to a formal 
Taylor series expansion in powers of (Za)'^ of the QED one-loop self-energy correction, which is expressed |8j: 



AE^ 



d^k i^D , 



'{2tt)^J fc2 + 



- (^^{r) |e*'^"-7^5j.(r,r',i?o-A:°)7^e-^''"-'| Vo^(r')) - Sm (^o'^lV'o'') 



(25) 



here 5i?(r,r',iJ) is the Dirac-Coulomb propagator and ^pP is the Dirac wave function. These extra terms do not 



appear in asymptotic expansion of the integrand in Eq. (25) and should be withdrawn. 



Thus our definition coincides with that of ^12 (see prescriptions after Eq. (3.16) in that Ref.^ 



III. NUMERICAL SCHEME 



Here we present a numerical scheme to evaluate the Bethe logarithm ([8|) and its relativistic corrections ( 24 ) for the 
two-center Coulomb problem. 



A. Variational expansion 

A variational expansion for the electronic wave function is taken in the form {Zi ^ Z2) 

00 



*„(r) = e""'^rl"l' 



(26) 



where r is a distance from the electron to the z-axis and 
1 



r = ^\/2^f '^i + 2^^?^?^ + 2r2i?2 -,,-,^ 



-m. 



For Zi = Z2 the variational wave function should be symmetrized 



— pi'm.vJ.™\ 



*(ri,r2) = e' -■-r'"' > o,;(e 



E^»(^ 



-OLiTx—liir^ 



±e 



-ftri— Qir2^ 



(27) 







R (in au) 



1,0 

R (in au) 



FIG. 4: The nonrelativistic Bethe logarithm for the ground (Iscr) electronic state for Z\ = Z2 = 1 (left) and Z\ = 2, Z2 = —1 
(right). Dashed lines show the vibrational wave functions for the ground state in HJ ion (left) and the (36, 34) state in ''He'^p 
(right). 

where (+) is used to get a gerade electronic state and (— ) is for an ungerade state, respectively. Parameters ai and 
/3i are generated in a quasi-random manner j23j 



«(» + l)\/p^ 



{A2-Ai)+Ai 



(28) 



Here [a;J designates the fractional part of x, Pa is a prime number, and [v4i, A2] is a real variational interval, which 
has to be optimized. Parameters Pi are obtained in a similar way. Details may be found elsewhere [551 HI] ■ 

All the integrands in Eqs. (l5|, ([l5]), (11a), (lib) and (13) have the form of a second-order perturbation expression, 
i.e. they involve an operator (Eq — H — k) . We thus diagonalize the matrix of the Hamiltonian for intermediate 
states to get a set of (pseudo)state wavefunctions ij^m and energies E„i. A similar approach was used to compute the 
nonrelativistic Bethe logarithm for the three-body Coulomb problem in [50] . 

The basis for intermediate states is constructed as follows. 

(1) We use a regular basis set (with regular values of the exponents a, (3), similar to that used for the initial state. 

(2) We build a special basis set with exponentially growing parameters for ri: 



,(0) 



Au 



,(0) 



^4")=r"Ai, 



A. 



4"^=r"A2, 



where r = A2/A1. Typically [^1,^2] = [2.5,4.5], and n„ 
interval fee [0,10"]. 



(29) 



5 — 7, which corresponds to the photon energy 



(3) We add a similar basis set for r2- Note that this last step may be omitted in the case Zi = Z2 where the basis is 
symmetrized. 



B. Nonrelativistic Bethe logarithm 



After expansion on the basis for intermediate state, the expression of the integrand becomes 



Pndik) - J2 



(V-olVI^™)' 



Eo — Em — k 
so the integral appearing in the low-energy part of the numerator is 

Eo — Em 



f kdkPUk) = Y.^^^\^\^- 
■Jo 



A-iEo- Em) In 



Eq — Em — A 



(30) 



(31) 





R (in au) 




R (in au) 



R (in au) 



FIG. 5: The four contributions to the relativistic Bathe logarithm for the ground (IscTg) electronic state, for Z\ = Z-i = 1. 
Dashed lines show their respective values for the Is state of a hydrogen atom [R] . 



It remains to calculate the matrix elements of the impulse operator. Its standard components are 



(32) 



Assuming, from now on, that "00 is a a-state, action on V'o of the impulse operator may be expressed as follows 



v[,'Vo 



R\ 1 



2 / n 



a., + z 



R\ 1 



2/ r2 



■00, 



V 






(33) 



Here z = (rj — r|)/(2i?). Using these relations then calculation of the matrix elements is straightforward 



C. Relativistic corrections to the Bethe logarithm 



We devote one paragraph to each term of Eq. ( 24 ) , giving additional details which are necessary for their numerical 

evaluation. Note that the terms /3}° and /3} are treated numerically in an independent way, which is a further reason 
for separating these two contributions. 



1. For /3}° , the integrand can be written as 



To evaluate (V'i|p^|V'i)i it is convenient to use the following identity (see Eq. (16l) 

p2Vi = 2[(So-V^-fc)V'i-V0o]. 



(34) 



(35) 





R (in au) 



FIG. 6: (Left) The relativistic Bethe logarithm for the ground (Isag) electronic state, for Zi — Z2 = 1. The dashed line is 
the jC{1s) for the Is state of a hydrogen atom [12]. (Right) The same data, but normalized by the delta function distribution; 
N{R) — Atv (Zi(5(ri) + Z25{r2)). The dashed line is the vibrational wave function for the ground state of H^. 



Then 

(V-ilP^l^i) = 4{ij,\{Eo-V-kf\^i) - A{tPi\{Eo-V~k)\ijo) + (V^) 
and, for arbitrary fc, 'ipi{k) may be expressed: 

I V'm)(V'm I V|^0) 



V'l(fc)-E 



Eo — Em — k 

,(lb)/ 



2. In order to get Pre [k) we first solve the equation 
It can be shown that ^b behaves at small ri (or r2) as 

z, zl 



'>pB{ri,r2) 



2 Inr-i j ^a{ri,r2) + ^^ 



Zo Zr, 



lnr2 'ipoiri,r2) + '0s(''i, ?'2) 



^4ri 2 J ""' "' \4r2 2 

where '>l)B{f 1,^2) is a regular function. Then equation for V'-b(^15^2) niay be written 



(^0 - H) i^i 



Hf 



HB)]i'0 + 



. Zl Z2 Zl Z2 

_H, In 7'i In To 

4ri 4r2 2 2 



V'o 



Thus substituting '0b(''i,''2) from Eq. (39 1 into Eq. (lib), one gets 
p(")(fc)^2((|L 



^' ' §-^-^lnri-^lnr2]QpiE,-H-k)-'p 



2UB\Qp{Eo~H-k) 'pUo 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 



The derivation of matrix elements involving logarithms may be found in the Appendix. 

3. For the evaluation of F„q(fc), Eq. (151 must be transformed to separate contributions from operators of different 
ranks. Denoting Ag = {Eq — H — k)^ , 



P,^g{k) ^-^ I d^n {6'^ - n'n^) n^n"" j k^dk {(pV A,rV) - (pV'r" A,pJ)} 



(42) 
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R 


^(.) 


^['' 


/32 


A 


£ 


C/N{R) 


0.2 


199.124 


-362.792 


-718.74 


500.9 


-381.5 


-9.4774 


0.3 


71.425 


-155.864 


-408.41 


244.09 


-248.76 


-7.7149 


0.4 


26.866 


-73.891 


-293.55 


153.814 


-186.76 


-7.0981 


0.5 


8.5247 


-36.4368 


-232.469 


111.218 


-149.163 


-6.8191 


0.6 


0.3419 


-17.8754 


-192.985 


86.9905 


-123.529 


-6.6763 


0.7 


-3.4173 


-8.19603 


-164.670 


71.3995 


-104.884 


-6.5980 


0.8 


-5.1044 


-2.99381 


-143.164 


60.5041 


-90.7581 


-6.5534 


0.9 


-5.7669 


-0.16341 


-126.255 


52.4484 


-79.7371 


-6.5273 


1.0 


-5.94144 


1.36042 


-112.653 


46.2540 


-70.9802 


-6.5145 


1.1 


-5.85242 


2.14699 


-101.528 


41.3527 


-63.8807 


-6.5084 


1.2 


-5.64166 


2.50981 


-92.3095 


37.3908 


-58.0506 


-6.5071 


1.3 


-5.37875 


2.62712 


-84.5938 


34.1346 


-53.2109 


-6.5097 


1.4 


-5.10067 


2.60269 


-78.0750 


31.4227 


-49.1503 


-6.5145 


1.5 


-4.82677 


2.49792 


-72.5270 


29.1400 


-45.7159 


-6.5212 


1.6 


-4.56684 


2.34961 


-67.7744 


27.2013 


-42.7903 


-6.5294 


1.7 


-4.32531 


2.17994 


-63.6799 


25.5425 


-40.2828 


-6.5387 


1.8 


-4.10371 


2.00222 


-60.1351 


24.1145 


-38.1221 


-6.5489 


1.9 


-3.90197 


1.82441 


-57.0531 


22.8786 


-36.2520 


-6.5598 


2.0 


-3.71921 


1.65108 


-54.3637 


21.8045 


-34.6274 


-6.5712 


2.1 


-3.55413 


1.48479 


-52.0099 


20.8676 


-33.2117 


-6.5831 


2.2 


-3.40532 


1.32679 


-49.9446 


20.0480 


-31.9752 


-6.5953 


2.3 


-3.27133 


1.17755 


-48.1290 


19.3294 


-30.8933 


-6.6077 


2.4 


-3.15078 


1.03708 


-46.5306 


18.6985 


-29.9458 


-6.6203 


2.5 


-3.04237 


0.90509 


-45.1223 


18.1441 


-29.1155 


-6.6329 


2.6 


-2.94493 


0.78117 


-43.8811 


17.6566 


-28.3883 


-6.6456 


2.7 


-2.85740 


0.66481 


-42.7875 


17.2283 


-27.7518 


-6.6583 


2.8 


-2.77881 


0.55550 


-41.8248 


16.8522 


-27.1958 


-6.6709 


2.9 


-2.70829 


0.45272 


-40.9786 


16.5527 


-26.6814 


-6.6758 


3.0 


-2.64510 


0.35594 


-40.2365 


16.2347 


-26.2909 


-6.6955 


3.2 


-2.53796 


0.17881 


-39.0229 


15.7664 


-25.6156 


-6.7192 


3.4 


-2.45271 


0.02095 


-38.1126 


15.4188 


-25.1256 


-6.7416 


3.6 


-2.38557 


-0.12013 


-37.4498 


15.1695 


-24.7860 


-6.7624 


3.8 


-2.33346 


-0.24634 


-36.9899 


15.0008 


-24.5689 


-6.7815 


4.0 


-2.29385 


-0.35913 


-36.6963 


14.8980 


-24.4513 


-6.7985 


4.2 


-2.26461 


-0.45961 


-36.5387 


14.8490 


-24.4139 


-6.8135 


4.4 


-2.24395 


-0.54861 


-36.4913 


14.8435 


-24.4404 


-6.8264 


4.6 


-2.23034 


-0.62682 


-36.5320 


14.8726 


-24.5165 


-6.8370 


4.8 


-2.22246 


-0.69487 


-36.6417 


14.9288 


-24.6303 


-6.8456 


5.0 


-2.21917 


-0.75336 


-36.8039 


15.0053 


-24.7712 


-6.8523 



TABLE I: Relativistic Bethe logarithm for the ground (1s(Tj) electronic state, for Zi — Z2 ~ 1. 



and one obtains 
3 



Stt 



d^n {S'^ - n'n^) n'ri"" / k^dk (pV' A, r"p^) 



= / k^dk 
^ I k^dk 



3 

10 
9 
20 



S 



(2) 



A,S, 



(2) 



p X r] As [r X p] 



^f]'^^f)-4([pxr]A, 



[rxp] 



(43) 



Stt 



dn^ {S'^ - n'n^) n^n"" / k^dk {p'r^r"' A,p^) 



= I k^dk \'^{p' r^Asp') - ^{{pry Asp^) \ . 



(44) 
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0,5 1,0 

R (in au) 



FIG. 7: (Left) The relativistic Bethe logarithm for the ground (Iscr) electronic state for Z\ =2 and Z2 ~ —1. (Right) The same 
data, but normalized by the delta function distribution; N{R) = 47r (^ZfS{ri) + Z2(5(r2)). The dashed line is the vibrational 
wave function for the (36, 34) state of *He"'"p. 



Here the quadrupole operator is defined as 



S. 



(2) _ 



1 



r,pj +rjp,- -(5y (r -p) 



and its standard components are 

Sq^^ = Si^J = 3 C^zp^ - xp^ - ypy) 

'^'±2 = ^ {S'i} " 4y^ ± '^iS'^xy) = —f^ [xPx -yPy±i {xPy + VPa;)] 

Matrix elements can be computed from the relations (valid for a a state): 



■-yz J =T^ [xpz +zpx±i [yp^ + zpy) 



S'o^ = 3 (2zp2+r+p_+r_p4 



2z 



^+f).^^- 



4^N;^(^...+^..)-±^--^((2.+ f);^a.. 



R\ 1 „ 

2z-f)^a.„ 



- r- I — dr, -\ 9^2 

'T 7-2 



;.(2) 
'±2 



V^ 



r±p± 



V6 V''i f2 



Here r± and p± are standard components of vector operators r and p, respectively. 
Finally, the matrix elements of the [r x p] operator can be obtained from 

- . ire±^^ R ( I ^ 1 ^ 



^(2) 



4. In the expression of Pre (k) in Eq. (13 1, the two terms can be respectively written as 



and 



2({Eo-V)p^Eo-H-k) >' 



U(z.^^^zJ^^]iE.-H-krp^ 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 
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R 


/?i 


132 


/?3 


C 


C/N{R) 


0.20 


29.65 


-584.1 


62.514 


-491.9 


-8.647 


0.25 


17.558 


-799.9 


159.373 


-622.9 


-9.1887 


0.30 


1.741 


-998.9 


236.611 


-760.5 


-9.5198 


0.35 


-15.627 


-1191.9 


305.525 


-902.0 


-9.7143 


0.40 


-33.158 


-1382.7 


369.822 


-1046.0 


-9.8402 


0.45 


-49.903 


-1566.8 


430.590 


-1186.1 


-9.8973 


0.50 


-65.195 


-1743.3 


487.854 


-1320.6 


-9.9213 


0.55 


-78.929 


-1909.9 


541.306 


-1447.5 


-9.9272 


0.60 


-91.000 


-2064.5 


590.636 


-1564.8 


-9.9209 


0.65 


-101.172 


-2207.3 


635.664 


-1672.8 


-9.9137 


0.70 


-109.370 


-2335.7 


676.345 


-1768.7 


-9.8946 


0.75 


-116.075 


-2451.0 


712.787 


-1854.2 


-9.8745 


0.80 


-121.547 


-2553.7 


745.209 


-1930.0 


-9.8546 


0.85 


-125.908 


-2645.8 


773.901 


-1997.8 


-9.8408 


0.90 


-127.846 


-2727.0 


799.189 


-2055.6 


-9.8199 


0.95 


-130.235 


-2797.5 


821.408 


-2106.3 


-9.8013 


1.00 


-131.770 


-2860.8 


840.902 


-2151.6 


-9.7898 


1.10 


-133.766 


-2963.8 


872.905 


-2224.6 


-9.7653 


1.20 


-134.563 


-3041.6 


897.365 


-2278.7 


-9.7402 


1.30 


-134.679 


-3101.3 


916.090 


-2319.8 


-9.7204 


1.40 


-134.430 


-3147.1 


930.468 


-2351.0 


-9.7042 


1.50 


-134.006 


-3182.5 


941.556 


-2374.9 


-9.6915 


1.60 


-133.517 


-3209.9 


950.148 


-2393.2 


-9.6810 


1.70 


-133.025 


-3231.3 


956.847 


-2407.4 


-9.6727 


1.80 


-132.561 


-3248.2 


962.102 


-2418.6 


-9.6664 


1.90 


-132.140 


-3261.4 


966.255 


-2427.2 


-9.6607 


2.00 


-131.766 


-3272.0 


969.562 


-2434.2 


-9.6565 



TABLE II: Relativistic Bethe logarithm for the ground (Isa) electronic state, for Zi — 2 and Z2 = —1. 

We now show that the second term does not contribute here. Since right-hand transition does not change spin, then 
on the left-hand side we may keep only those terms, which contain a^'- 

TU^^'-^(^ + ^)iEo-H-kr^A (51) 



This contribution results in fine splitting of the main line, which is spin-dependent, and thus will not be considered 
here. 

IV. RESULTS 

A. Nonrelativistic Bethe logarithm 

First we perform calculations of the nonrelativistic Bethe logarithm for the two-center Coulomb problem for two 
cases: Zi = Z2 = 1 and Zi = 2, Z2 = —1. Results are plotted in Fig. [4] The numerical data for Pnr{R) is as accurate 
as 10 significant digits and has been obtained for a large range of internuclear bond lengths i? = — 7 a.u. with a 
step of 0.05 a.u. For the hydrogen molecular ion our results are in a good agreement with previous calculations |25j . 
The complete Tables are too lengthy to be reported here and may be found in the Supplemental Material [26] . 

B. Relativistic corrections to the Bethe logarithm 

The numerical scheme has been tested on the ground state of the hydrogen atom. Our results are 

/3i(ls) = -3.26821319, 
^2(ls) = -40.647026 693. 
^3(ls) = 16.655 330 436, 
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They are in perfect agreement with the ones of |12) . 

The main results of this work are presented in Tables IT] and lllj Particular behaviour of different components of the 
relativistic Bethe logarithm £, namely P^ \ Pi , 1^2, P3, for a case of a hydrogen molecular ion are shown on Fig. ^ 
On Figures Ig] and [7] our final results for the low energy contribution C{R) are plotted. 

The numerical calculations have been performed in multiple precision arithmetic with 64 decimal digits. Special 
care has been required for the nonrelativistic quadrupole contribution at small R, in this case a 96 decimal digit 
arithmetic has been used. 

It is worth noting that in the region where the wave function is essentially nonzero the quantity C/N{R) is about 
constant: C/N{R) sa 6.5 for Zi = Z2 = 1 and C/N{R) w 9.7 for Zi == 2, Z2 = -1. That may help in a qualitative 
estimate within 2-3 digit accuracy of the one-loop self-energy contribution to the transition frequencies of not only 
hydrogen molecular ions but also of the neutral hydrogen molecule H2 . 

In conclusion, we have computed the low-energy part of the ma'^ order self-energy correction for the two center 
problem, with a numerical accuracy that exceeds 3 significant digits, in the whole range of internuclear distances 
R G [0.2,7]. Calculation of the high energy part is under consideration now. A complete result then should be 
averaged over vibrational wave functions to get proper correction for ro- vibrational transition frequencies at the ma^ 
order. Taking into account as well the vacuum polarization and two-loop electron self-energy contributions at this 
order we expect that a relative precision for transition frequencies will be at a level of 10"^" or better both for 
hydrogen molecular ions and for antiprotonic helium. 
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Appendix: Two-center matrix elements involving Inri or lnr2 

In this Appendix we want to consider an analytical evaluation of the two-center integrals of type [22] , which contain 
either Inri or lnr2 as a multiplier. We start from the basic integral: 



-ari —/3r2 



rinri,oo(a,/3) = / dridr2 Inri • e 

Using identities 

1"°° 1 

/ drlnr e"'^'' = - [Ei{Rj) +lnR- e"'^^] 



(52) 



and 



/ drlnr e''"' = [Ei(i?7) + In 7 + In R ■ e-''^ + 7^] 



one may get 



/■ri+R 

dn Inrie"""^! / drs e"'''^^ 

R Jri-R 



\nrie-°"'' ( 



art l^-Piri-R) _ g-/3(ri+fl) 



(53a) 



O0R „ p-PR 



Ei(i?(a-f/3)) +\nR- 
a + P 



,-ia+l3)R 
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and 



R pR+ri 

(iri Inri e^"'^^ / dr2 e 

JR-ri 

R 



'0r2 _ 



1 

p-PR 



dri 



-an l^-P(R-ri) _g-/3(i?+ri) 



-PR 



/3 



Inrie-"''! f 
R (.it: 

Ei(i?(a+/?)) + In (a+/3) + In i? • e-("+'^)'^ + 7^ 
a + /3 

Ei(i?(a-/3)) + In {a~ {3) + Ini? • g-^""'')^ + 7e 



(53b) 



a — /3 
Summing up the two contributions we obtain a final expression 

47r e-/^-" 7b + e^"^ In i? 27r e^^-^ Ei(i?(a+/3)) + e-^^-" In (a + /3) 



Tin n, 00 (0:7/^) 



i? 



a2-/3 



RP a + 13 

27^ e-'^^ [Ei(fl(a-^)) +ln(a-/3)] 
'liJ3 a- 13 



(54) 



To generate other integrals one may use 
da) \ d(3 



^Inn, kl{(X,0} 



rinri,Oo(a,/9). 



(55) 
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In this Supplemental Material, we give the nonrclativistic Bethe logarithm (Table I) and more extensive results for 
its relativistic corrections (Tables II and III), both for the ground electronic state of H^ and Hc^p. These additional 
data arc helpful for computational use to perform more accurately averaging over vibrational wave functions. 

Nonrelativistic Bethe logarithm 



TABLE III: The nonrelativistic Bethe logarithm, /3„ 



for HJ and He p. 



R (in a.u.) 


^1 = 1,^2 = 1 


Zi^2, Z2 = -l 


0.00 


4.370422917 


2.984128556 


0.05 


3.960959074 


4.429093491 


0.10 


3.763208049 


4.695144972 


0.15 


3.627181789 


4.737059502 


0.20 


3.525244876 


4.712352798 


0.25 


3.445386029 


4.668731134 


0.30 


3.381060204 


4.622934909 


0.35 


3.328241427 


4.580920363 


0.40 


3.284256226 


4.544480350 


0.45 


3.247231763 


4.513745134 


0.50 


3.215803070 


4.488209081 


0.55 


3.188944303 


4.467169500 


0.60 


3.165865517 


4.449913692 


0.65 


3.145946182 


4.435793560 


0.70 


3.128690731 


4.424249332 


0.75 


3.113697833 


4.414810709 


0.80 


3.100638579 


4.407088769 


0.85 


3.089240633 


4.400764849 


0.90 


3.079276477 


4.395579241 


0.95 


3.070554542 


4.391320901 


1.00 


3.062912414 


4.387818567 


1.05 


3.056211557 


4.384933303 


1.10 


3.050333156 


4.382552371 


1.15 


3.045174825 


4.380584212 


1.20 


3.040647951 


4.378954395 


1.25 


3.036675548 


4.377602338 


1.30 


3.033190500 


4.376478663 


1.35 


3.030134120 


4.375543070 


1.40 


3.027454949 


4.374762617 


1.45 


3.025107756 


4.374110337 


1.50 


3.023052703 


4.373564122 


1.55 


3.021254630 


4.373105821 


1.60 


3.019682462 


4.372720508 


1.65 


3.018308688 


4.372395891 


1.70 


3.017108931 


4.372121836 


1.75 


3.016061565 


4.371889968 


1.80 


3.015147393 


4.371693363 


1.85 


3.014349362 


4.371526283 


1.90 


3.013652319 


4.371383970 


1.95 


3.013042797 


4.371262469 



TABLE III: (continued) 
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R (in a.u.) 


^1 = 1,^2 = 1 


Zi = 2, Z2 = -l 


2.00 


3.012508830 


4.371158489 


2.05 


3.012039785 


4.371069288 


2.10 


3.011626220 


4.370992576 


2.15 


3.011259757 


4.370926439 


2.20 


3.010932968 


4.370869276 


2.25 


3.010639275 


4.370819742 


2.30 


3.010372862 


4.370776708 


2.35 


3.010128599 


4.370739223 


2.40 


3.009901965 


4.370706488 


2.45 


3.009688991 


4.370677824 


2.50 


3.009486203 


4.370652660 


2.55 


3.009290569 


4.370630511 


2.60 


3.009099455 


4.370610966 


2.65 


3.008910587 


4.370593673 


2.70 


3.008722010 


4.370578334 


2.75 


3.008532059 


4.370564694 


2.80 


3.008339327 


4.370552535 


2.85 


3.008142638 


4.370541670 


2.90 


3.007941025 


4.370531937 


2.95 


3.007733705 


4.370523199 


3.00 


3.007520064 


4.370515335 


3.05 


3.007299631 


4.370508243 


3.10 


3.007072071 


4.370501832 


3.15 


3.006837162 


4.370496026 


3.20 


3.006594789 


4.370490756 


3.25 


3.006344924 


4.370485963 


3.30 


3.006087623 


4.370481595 


3.35 


3.005823009 


4.370477609 


3.40 


3.005551267 


4.370473962 


3.45 


3.005272634 


4.370470622 


3.50 


3.004987392 


4.370467557 


3.55 


3.004695862 


4.370464739 


3.60 


3.004398396 


4.370462145 


3.65 


3.004095374 


4.370459753 


3.70 


3.003787194 


4.370457545 


3.75 


3.003474273 


4.370455503 


3.80 


3.003157041 


4.370453612 


3.85 


3.002835935 


4.370451859 


3.90 


3.002511398 


4.370450231 


3.95 


3.002183876 


4.370448719 


4.00 


3.001853814 


4.370447311 


4.05 


3.001521655 


4.370445999 


4.10 


3.001187836 


4.370444776 


4.15 


3.000852787 


4.370443634 


4.20 


3.000516930 


4.370442567 


4.25 


3.000180678 


4.370441568 


4.30 


2.999844431 


4.370440633 


4.35 


2.999508578 


4.370439757 


4.40 


2.999173494 


4.370438934 


4.45 


2.998839540 


4.370438162 


4.50 


2.998507062 


4.370437437 


4.55 


2.998176392 


4.370436754 


4.60 


2.997847844 


4.370436112 


4.65 


2.997521719 


4.370435506 


4.70 


2.997198299 


4.370434936 


4.75 


2.996877852 


4.370434397 


4.80 


2.996560629 


4.370433889 


4.85 


2.996246864 


4.370433409 


4.90 


2.995936776 


4.370432955 


4.95 


2.995630567 


4.370432525 
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TABLE III: (continued) 



R (in a.u.) 


^1 = 1,^2 = 1 


Zi = 2, Z2 = -l 


5.00 


2.995328425 


4.370432118 


5.05 


2.995030520 


4.370431733 


5.10 


2.994737010 


4.370431368 


5.15 


2.994448035 


4.370431021 


5.20 


2.994163724 


4.370430693 


5.25 


2.993884189 


4.370430380 


5.30 


2.993609531 


4.370430084 


5.35 


2.993339837 


4.370429802 


5.40 


2.993075180 


4.370429533 


5.45 


2.992815623 


4.370429278 


5.50 


2.992561218 


4.370429035 


5.55 


2.992312003 


4.370428803 


5.60 


2.992068009 


4.370428583 


5.65 


2.991829253 


4.370428372 


5.70 


2.991595748 


4.370428171 


5.75 


2.991367492 


4.370427980 


5.80 


2.991144480 


4.370427797 


5.85 


2.990926695 


4.370427622 


5.90 


2.990714115 


4.370427455 


5.95 


2.990506709 


4.370427295 


6.00 


2.990304442 


4.370427142 


6.05 


2.990107272 


4.370426996 


6.10 


2.989915150 


4.370426856 


6.15 


2.989728023 


4.370426722 


6.20 


2.989545835 


4.370426593 


6.25 


2.989368523 


4.370426470 


6.30 


2.989196022 


4.370426352 


6.35 


2.989028262 


4.370426239 


6.40 


2.988865171 


4.370426130 


6.45 


2.988706674 


4.370426025 


6.50 


2.988552692 


4.370425925 


6.55 


2.988403148 


4.370425829 


6.60 


2.988257957 


4.370425736 


6.65 


2.988117039 


4.370425647 


6.70 


2.987980307 


4.370425562 


6.75 


2.987847676 


4.370425479 


6.80 


2.987719061 


4.370425400 


6.85 


2.987594373 


4.370425324 


6.90 


2.987473527 


4.370425251 


6.95 


2.987356433 


4.370425180 


7.00 


2.987243004 


4.370425112 



Relativistic correction to the Bethe logarithm 



TABLE IV: Relativistic Bethe logarithm for the ground (Isag) electronic 
state for Zi = Z2 = 1- 



R 


^(.) 


Pf' 


/32 


/?3 


C 


C/N{R) 


0.2 


199.124 


-362.792 


-718.74 


500.9 


-381.5 


-9.4774 


0.3 


71.425 


-155.864 


-408.41 


244.09 


-248.76 


-7.7149 


0.4 


26.866 


-73.891 


-293.55 


153.814 


-186.76 


-7.0981 


0.5 


8.5247 


-36.4368 


-232.469 


111.218 


-149.163 


-6.8191 


0.6 


0.3419 


-17.8754 


-192.985 


86.9905 


-123.529 


-6.6763 


0.7 


-3.4173 


-8.19603 


-164.670 


71.3995 


-104.884 


-6.5980 


0.8 


-5.1044 


-2.99381 


-143.164 


60.5041 


-90.7581 


-6.5534 


0.9 


-5.7669 


-0.16341 


-126.255 


52.4484 


-79.7371 


-6.5273 


1.0 


-5.94144 


1.36042 


-112.653 


46.2540 


-70.9802 


-6.5145 


1.1 


-5.85242 


2.14699 


-101.528 


41.3527 


-63.8807 


-6.5084 



TABLE IV: (continued) 
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R 


^(.) 


/3f' 


/32 


133 


C 


C/N{R) 


1.2 


-5.64166 


2.50981 


-92.3095 


37.3908 


-58.0506 


-6.5071 


1.3 


-5.37875 


2.62712 


-84.5938 


34.1346 


-53.2109 


-6.5097 


1.4 


-5.10067 


2.60269 


-78.0750 


31.4227 


-49.1503 


-6.5145 


1.5 


-4.82677 


2.49792 


-72.5270 


29.1400 


-45.7159 


-6.5212 


1.6 


-4.56684 


2.34961 


-67.7744 


27.2013 


-42.7903 


-6.5294 


1.7 


-4.32531 


2.17994 


-63.6799 


25.5425 


-40.2828 


-6.5387 


1.8 


-4.10371 


2.00222 


-60.1351 


24.1145 


-38.1221 


-6.5489 


1.9 


-3.90197 


1.82441 


-57.0531 


22.8786 


-36.2520 


-6.5598 


2.0 


-3.71921 


1.65108 


-54.3637 


21.8045 


-34.6274 


-6.5712 


2.1 


-3.55413 


1.48479 


-52.0099 


20.8676 


-33.2117 


-6.5831 


2.2 


-3.40532 


1.32679 


-49.9446 


20.0480 


-31.9752 


-6.5953 


2.3 


-3.27133 


1.17755 


-48.1290 


19.3294 


-30.8933 


-6.6077 


2.4 


-3.15078 


1.03708 


-46.5306 


18.6985 


-29.9458 


-6.6203 


2.5 


-3.04237 


0.90509 


-45.1223 


18.1441 


-29.1155 


-6.6329 


2.6 


-2.94493 


0.78117 


-43.8811 


17.6566 


-28.3883 


-6.6456 


2.7 


-2.85740 


0.66481 


-42.7875 


17.2283 


-27.7518 


-6.6583 


2.8 


-2.77881 


0.55550 


-41.8248 


16.8522 


-27.1958 


-6.6709 


2.9 


-2.70829 


0.45272 


-40.9786 


16.5527 


-26.6814 


-6.6758 


3.0 


-2.64509 


0.35594 


-40.2365 


16.2347 


-26.2909 


-6.6955 


3.1 


-2.58852 


0.26478 


-39.5878 


15.9839 


-25.9276 


-6.7075 


3.2 


-2.53795 


0.17881 


-39.0229 


15.7664 


-25.6156 


-6.7192 


3.3 


-2.49285 


0.09764 


-38.5337 


15.5791 


-25.3498 


-6.7306 


3.4 


-2.45270 


0.02094 


-38.1126 


15.4188 


-25.1256 


-6.7416 


3.5 


-2.41708 


-0.05155 


-37.7533 


15.2831 


-24.9389 


-6.7522 


3.6 


-2.38556 


-0.12013 


-37.4498 


15.1695 


-24.7860 


-6.7624 


3.7 


-2.35780 


-0.18499 


-37.1970 


15.0761 


-24.6637 


-6.7722 


3.8 


-2.33346 


-0.24634 


-36.9899 


15.0008 


-24.5689 


-6.7815 


3.9 


-2.31223 


-0.30433 


-36.8243 


14.9419 


-24.4990 


-6.7903 


4.0 


-2.29385 


-0.35913 


-36.6963 


14.8980 


-24.4513 


-6.7985 


4.1 


-2.27805 


-0.41085 


-36.6028 


14.8674 


-24.4243 


-6.8065 


4.2 


-2.26461 


-0.45960 


-36.5387 


14.8490 


-24.4139 


-6.8135 


4.3 


-2.25331 


-0.50549 


-36.5027 


14.8414 


-24.4201 


-6.8202 


4.4 


-2.24395 


-0.54860 


-36.4913 


14.8435 


-24.4404 


-6.8264 


4.5 


-2.23635 


-0.58902 


-36.5019 


14.8542 


-24.4730 


-6.8320 


4.6 


-2.23034 


-0.62682 


-36.5320 


14.8726 


-24.5165 


-6.8370 


4.7 


-2.22576 


-0.66208 


-36.5793 


14.8978 


-24.5694 


-6.8416 


4.8 


-2.22246 


-0.69487 


-36.6417 


14.9288 


-24.6303 


-6.8456 


4.9 


-2.22031 


-0.72527 


-36.7172 


14.9649 


-24.6979 


-6.8492 


5.0 


-2.21917 


-0.75335 


-36.8039 


15.0053 


-24.7712 


-6.8523 


5.1 


-2.21893 


-0.77921 


-36.9002 


15.0493 


-24.8490 


-6.8549 


5.2 


-2.21949 


-0.80291 


-37.0044 


15.0964 


-24.9304 


-6.8571 


5.3 


-2.22073 


-0.82454 


-37.1151 


15.1459 


-25.0145 


-6.8588 


5.4 


-2.22257 


-0.84421 


-37.2308 


15.1972 


-25.1004 


-6.8602 


5.5 


-2.22492 


-0.86200 


-37.3505 


15.2499 


-25.1875 


-6.8612 


5.6 


-2.22770 


-0.87800 


-37.4729 


15.3036 


-25.2750 


-6.8619 


5.7 


-2.23084 


-0.89231 


-37.5970 


15.3578 


-25.3624 


-6.8623 


5.8 


-2.23427 


-0.90504 


-37.7220 


15.4121 


-25.4493 


-6.8624 


5.9 


-2.23794 


-0.91629 


-37.8470 


15.4662 


-25.5350 


-6.8623 


6.0 


-2.24178 


-0.92615 


-37.9714 


15.5199 


-25.6194 


-6.8619 


6.1 


-2.24576 


-0.93472 


-38.0943 


15.5729 


-25.7019 


-6.8613 


6.2 


-2.24982 


-0.94211 


-38.2154 


15.6250 


-25.7824 


-6.8606 


6.3 


-2.25394 


-0.94839 


-38.3342 


15.6760 


-25.8606 


-6.8597 


6.4 


-2.25808 


-0.95368 


-38.4501 


15.7256 


-25.9363 


-6.8587 


6.5 


-2.26220 


-0.95805 


-38.5630 


15.7739 


-26.0093 


-6.8576 


6.6 


-2.26629 


-0.96158 


-38.6725 


15.8207 


-26.0797 


-6.8563 


6.7 


-2.27032 


-0.96439 


-38.7784 


15.8658 


-26.1473 


-6.8551 


6.8 


-2.27427 


-0.96652 


-38.8806 


15.9094 


-26.2120 


-6.8537 


6.9 


-2.27814 


-0.96805 


-38.9789 


15.9512 


-26.2739 


-6.8523 


7.0 


-2.28190 


-0.96905 


-39.0732 


15.9913 


-26.3328 


-6.8509 
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TABLE V: Relativistic Bethe logarithm for the ground (Isa) electronic 
state, for Zi = 2 and Z2 = —1. 



R 


/?i 


132 


/33 


C 


C/N{R) 


0.10 


192.1 


-102. 


-381.08 


-290.9 


-7.3871 


0.15 


38.01 


-288. 


-81.305 


-331.2 


-6.9913 


0.20 


29.65 


-584.1 


62.514 


-491.9 


-8.6471 


0.25 


17.558 


-799.9 


159.373 


-622.9 


-9.1887 


0.30 


1.741 


-998.9 


236.611 


-760.5 


-9.5198 


0.35 


-15.627 


-1191.9 


305.525 


-902.0 


-9.7143 


0.40 


-33.158 


-1382.7 


369.822 


-1046.0 


-9.8402 


0.45 


-49.903 


-1566.8 


430.590 


-1186.1 


-9.8973 


0.50 


-65.195 


-1743.3 


487.854 


-1320.6 


-9.9213 


0.55 


-78.929 


-1909.9 


541.306 


-1447.5 


-9.9272 


0.60 


-91.000 


-2064.5 


590.636 


-1564.8 


-9.9209 


0.65 


-101.172 


-2207.3 


635.664 


-1672.8 


-9.9137 


0.70 


-109.370 


-2335.7 


676.345 


-1768.7 


-9.8946 


0.75 


-116.075 


-2451.0 


712.787 


-1854.2 


-9.8745 


0.80 


-121.547 


-2553.7 


745.209 


-1930.0 


-9.8546 


0.85 


-125.908 


-2645.8 


773.901 


-1997.8 


-9.8408 


0.90 


-127.846 


-2727.0 


799.189 


-2055.6 


-9.8199 


0.95 


-130.235 


-2797.5 


821.408 


-2106.3 


-9.8013 


1.00 


-131.770 


-2860.8 


840.902 


-2151.6 


-9.7898 


1.05 


-132.958 


-2916.3 


857.985 


-2191.2 


-9.7795 


1.10 


-133.766 


-2963.8 


872.905 


-2224.6 


-9.7653 


1.15 


-134.277 


-3005.3 


885.953 


-2253.6 


-9.7521 


1.20 


-134.563 


-3041.6 


897.365 


-2278.7 


-9.7402 


1.25 


-134.682 


-3073.4 


907.349 


-2300.7 


-9.7297 


1.30 


-134.679 


-3101.3 


916.090 


-2319.8 


-9.7204 


1.35 


-134.586 


-3125.7 


923.750 


-2336.5 


-9.7119 


1.40 


-134.430 


-3147.1 


930.468 


-2351.0 


-9.7042 


1.45 


-134.232 


-3165.9 


936.368 


-2363.7 


-9.6974 


1.50 


-134.006 


-3182.5 


941.556 


-2374.9 


-9.6915 


1.55 


-133.765 


-3197.1 


946.122 


-2384.7 


-9.6861 


1.60 


-133.517 


-3209.9 


950.148 


-2393.2 


-9.6810 


1.65 


-133.268 


-3221.3 


953.703 


-2400.8 


-9.6768 


1.70 


-133.025 


-3231.3 


956.847 


-2407.4 


-9.6727 


1.75 


-132.788 


-3240.3 


959.631 


-2413.4 


-9.6696 


1.80 


-132.561 


-3248.2 


962.102 


-2418.6 


-9.6664 


1.85 


-132.344 


-3255.2 


964.299 


-2423.2 


-9.6635 


1.90 


-132.140 


-3261.4 


966.255 


-2427.2 


-9.6607 


1.95 


-131.947 


-3267.0 


968.001 


-2430.9 


-9.6585 


2.00 


-131.766 


-3272.0 


969.562 


-2434.2 


-9.6565 



